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ABSTRACT 


An examination of the potential flow computer code VSAERO to model 
leading edge separation over a delta wing. Recent improvements to the code 
suggest that it may be capable of predicting pressure coefficients on the 
body. Investigation showed that although that code does predict the vortex 
roll-up, the pressure coefficients have significant error. The program is 
currently unsatisfactory, but with some additional development it may 
become a useful tool for this application. 
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NOMENCLATURE 


Local span 
Chord length 

Differential surface element 

Unit normal vector 

Arbitrary point in space 

Vector between P and surface element dS 

Surface of the arbitrary body 

Imaginary surface at infinity 

Wake surface 

Velocity vector 

Chordwise coordinate 

Spanwise coordinate 

Coordinate perpendicular to x-y plane 

Velocity potential 

Doublet singularity strength per unit area 

Source singularity strength per unit area 

subscripts 

Interior region 

Lower wake surface 

Refers to point P 

Upper wake surface 

Free stream conditions 
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INTRODUCTION 


Panel methods have been primarily used for the determination of a 
flow field about a two or three dimensional body. The body is represente as 
being made up of a set of panels with a distribution of singularities across 
the set. If the singularities are constant within a panel a low order code is 
the result. If the singularities are variable within a panel a high order code 
is the result. It has been shown that low order panel codes can achieve good 
results if the panel density is sufficiently high at a substantial savings in 
computing time. The singularities can be sources, sinks, doublets, vortices, 
and any combination. However, if the body is to generate lift the panels 
must be either doublets or vortices. 

The strength of panel methods is the capacity to model complex, 
arbitrary configurations. There has been some interest in using the low 
order panel code VSAERO to model a highly swept delta wing at high 
angles of attack. If this can be achieved it will be very useful in the design 
and analysis of high performance aircraft and extra-atmospheric vehicles 
such as the National Aerospace Plane. 

When the leading edge of a delta wing is sharp, the flow separates 
from the wing, forms a spiral vortex sheet and then reattaches to the wing. 
The vortex sheet induces increased velocities on the surface of the wing. 

As a result the lift coefficient is higher than for attached flow for the same 
flight conditions. This incremental lift is called the vortex or non-linear 
lift and is highly desirable. 

The difficulty to the problem is that panel codes are based on the 
invisdd flow equations, and model the effects of viscosity by the means of a 
wake composed of source panels. Potential flow codes historically have 
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not produced good results when separation has a significant effect. 
However, preliminary studies have shown some promise. 
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THEORY 


VSAERO uses the low order potential flow equations developed by 
Lamb (ref. 1). As we know from potentioal flow theory, any flow solution 
must satisfy the Lapace's Equation. Lamb assumes that a flow would have 
two solutions <j> and <j)' such that: 

V $ = V <j>' = 0 (1) 


The Divergence theorem states that 


J V 2 <t>dV = -/ s V<|>dS 


( 2 ) 


applying the two equations, and multiplying by <J>’ results in following 
equation: 



V 2 <j>dV = / <J)' V<{>dS = 


0 . 


( 3 ) 


<t> and 4»' can be reversed to yield 

/♦•V»dS = /*VfdS.O. 

o b (4) 

Lamb chooses <f>’ to be 1/r which satisfies the Laplacian and is a relatively 
simple function and the equation becomes 
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( 5 ) 




dS. 


However, a singularity occurs at r = 0. By describing a small sphere with 
an incremental area of de a about point P the singularity is excluded and 
the equation becomes: 



If we evaluate <J> at point P the equation becomes: 




.-4 

4iy 


<l>v(i|dS + 


-f 1 

4n J r 


V<j>dS. 


(7) 


This is 4> at any point P of the fluid in terms of <j> and V<}> at the boundary. 
The first integral is the disturbance potential of a distribution of double 
sources or doublets with axes normal to the surface having a density <}> per 
unit area, while the second integral is the disturbance potential of a 
distribution of simple sources having a density of V<}> per unit area. If P is 
external to the surface then: 




0 



<t>iV(l)dS 

\'l 



(8) 
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Subtracting equation 8 from equation 7 results in an equation for the 
velocity potential of the fluid both internal and external to the surface. As 
shown in Figure 1 the internal flow is the flow of interest, and the external 
flow is fictitious. The equation becomes: 


<t>D = 


— I (<t> - <t>i)n*v(!)dS + — I inM - V<t>j)dS 
4nj s+ WfS. 4 ' r ' 4nJ s+ms J 1 * 


(9) 


VSAERO assumes that S~ can be considered a large sphere centered 
on P such that: 



Since the disturbance potential at infinity is zero, the velocity potential at 
infinity must that due to the onset flow. For W the upper and lower surfaces 
are infinitesimally close, and that the corresponding upper and lower 
elements can be combined <|>i = 0. In addition entrainment is ignored 
- 0 . 

These assumptions result in: 

♦Jw*. ■ / w (V - < n ^ • v(i) dw ♦ ^ 


The surface S does not have any simplifying assumptions, however, if the 
point P lies on S the integral becomes singular. P is excluded from the 


5 


surface S by describing a local hemispherical defromation, and <)> due to the 
local hemispherical deformation is: 



-Lj f (<t>-f)n .v(l)R 2 sinydrde 

■iff* 


1 (v<J) - V<|)’) R sinydydG. 


(ID 


The integrals are evaluated between zero and pi, and the limit as R 
approaches zero results in: 



( 12 ) 


If equations 9, 10, and 12 are combined appropriately we arrive at the 
equation for <j> at any point P based on the values of <}> on the surface, the 
equation is: 



The boundary condition used to solve for <J) P is the internal Dirichlet 
boundary condition. The total potential <j) is taken as being composed of the 
onset potential <)>« , and a disturbance potential <J>d = <J> - <j>«» . In order to 
minimize the size of the surface singularities the potential of the fictitious 
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flow is set equal to that of the onset flow. Using these boundary conditions 
and looking at the ficticitious flow equation 13 becomes: 


0 = 


_Lj^n.v(l) dS + _lJ s !„^-v^s 


(14) 


The doublet strength is defined as: 

4IT|i = <t) d = 4> -<!>*, 

The source strength is defined as: 

4IIa = -n • (v<J> - V<(>^ = V N - n • , 

and can be solved directly for c. For a more general case V N may have two 
components representing a boundary layer displacement and 
inflow/outflow. Knowing the source strengths, VSAERO then solves 
equation 14 for the unknown doublet strengths. 



If P is off the surface, then K is zero. If P is on an inside smooth suface, 
then K is -2k. If P is on an outside smooth suface, then K is 2k. If P is on a 
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crease , then K is equal to the solid angle. If the potential field has been 
computed for a grid, the velocity field can be calculated using a local 

differentiation scheme since: 

V P = -V<|) p , 

or calculated directly which is the current VSAERO methodwith the 
following equation: 
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DISCUSSION 


Initially the 3000 panel version of VSAERO was used in the attempt to 
model the leading edge separation problem. The literature suggested that 
the code would converge on a solution. However, it appears that the user 
would have to enter a wake shape very close to the actually shape. There is 
a preprocessor commercially available called VORSEP which is suppose to 
generate a wake close to actual that allows the code to converge to the 
solution. This preprocessor was not immediately available, and it was 
decided to examine a modified version of VSAERO that was thought to be 
capable modeling the problem. 

The modified version is time-dependent, and the wake is not 
specified. At time t = 0 the body does not have a wake. At time t = ti it has a 
wake one panel in length. The generation of the wake is shown in figure 23. 
The wake is generated by following the path of a particle originating on the 
separation line where the velocity of the particle at point P is equal to the 
velocity of the fluid at point P. In this manner it was expected that wake 
generated would follow the path of the vortex sheet, and produce reasonable 
pressure distributions. 

Be aware that the program has many variables which can be 
adjusted to modify the calculations. After an exhaustive parametric study 
it was determined that the program was unable to produce results to the 
desired accuracy. The wake shape generated was as far as the eye could 
tell correct. The pressure coefficients were incorrect in magnitude, but the 
shape of the curve of Cp versus span was correct. It was thought that the 
method of calculating the particle trajectory was too simple. Several 
different numerical schemes were initiated including a Runge-Kutte and 
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Adam's Predictor-Corrector methods. For each modification the same 
parametric study was performed. The results being a substantial increased 
in computing time without a significant improvement to the solution. 

The velocity of the particle has three components: the onset velocity, 
the velocity due to body, and the velocity due to wake. Of these three only the 
velocity calculations due to the wake was new coding. The lines making up 
a wake panel are considered to be vortex filaments. Since the velocity due to 
an inviscid, vortex filament approaches infinity as the distance from that 
filament, R, approaches zero, there is a cut off distance called RCUT. When 
R is less than RCUT the induced velocity is set to zero. This was modified to 
such that when R was less than RCUT the induced velocity became a linear 
function of R. In either version RCUT is a very important parameter. This 
modification had little effect on the results. 

At this it was decided to test the capability of the program by 
manually forcing the wake into the proper shape. This produced improved 
results the best results thus far. The results, shown in figures 2 through 8, 
are not satisfactory. The program cannot produce better results without 
substantial modification. 
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CONCLUSIONS 


VSAERO as it currently stands is inadequate for the modeling of leading 
edge separation. In addition to the poor correspondence with the Hummel 
wing, there is another problem. VSAERO assumes that the wake connects 
the body to a large spherical surface at infinity. Due to the neccessity of 
taking very small time increments in modeling of leading edge separation, 
the length of the wake is on the order of the maximum body dimension. 
However, this does not appear to be as important as the near field problem. 

The near field calculation coding is known to have problems, and in 
the investigationit was proven so. In order to obtain a reasonable flowfield 
about the body the near field calculations had to be bypass. Since the wake 
is in the near field, this has the potential for considerable error. If the near 
field problem can be corrected, VSAERO could be a useful tool in the 
analysis of leading edge separation in complex configurations. 
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Figure 2 Non-linear lift 
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Figure 3 Pressure coefficient vs. span 
at X=.3C 


VSAERO 


6 




y/b 


Hummel 

VSAERO 


Figure 4 


Pressure Coefficient vs. Span at 





Hummel 

VSAERO 


Figure 5 Pressure Coefficient vs. span at 
X=.7C 
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Figure 6 Pressure coefficient vs. span at 
X=.9C 
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